Interpreting scattering wave functions in the presence of energy-dependent 

interactions 
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In scattering theory, the squared relative wave function |^(q, r)| 2 is often interpreted as a weight, 
due to final-state interactions, describing the probability enhancement for emission with asymptotic 
relative momentum q. An equivalence relation also links the integral of the squared wave function 
over all coordinate space to the density of states. This relation, which plays an important role in 
understanding two-particle correlation phenomenology, is altered for the case where the potential is 
energy dependent, as is assumed in various forms of reaction theory. Here, the modification to the 
equivalence relation is derived, and it is shown that the squared wave function should be augmented 
by a additional factor if it is to represent the emission enhancement for final-state interactions. 
Examples with relativistic vector interactions, e.g., the Coulomb interaction, are presented. 



In many forms of reaction theory the square of the outgoing wave function, |0(q, r)| 2 , plays the role of a 
weight, enhancing the probability for emission from coordinate r into an asymptotic momentum state q. One 
on : such example is the Koonin equation [H, 0], used for two-particle correlations, 

C(q)= /d 3 rS(r)|<Xq,r)| 2 , (1) 



where q is the relative momentum and S(r) is the normalized probability for emitting two particles of the same 



momentum so that in their center-of-mass frame they are separated by r. Thus, the squared wave function 
describes the additional probability for emitting particles due to their final-state interactions. A second method 
for calculating the probability enhancement is to consider the correction to the relative density of states, which 
O^l ' is related to the phase shifts by the Beth-Uhlenbeck formula [!, 0, Hj] , 

m ; 

where Adn/de is the change of the density of states induced by the potential. The Beth-Uhlenbeck formula is 
I ■ only applicable for volumes V which are much larger than the range of the interaction or the inverse momentum 
f**"* \ 1/q. From this perspective, final-state interactions provide an extra weight for emission through the phase shift, 



where the denominator is the density of free states and q is the relative momentum. Since the volume V is 
in the denominator, this relation can be used to experimentally infer the size of the system. For the case of 
pion-proton correlations, a bump ensues in the correlation function for values of q corresponding to the invariant 
mass of the delta resonance. In this neighborhood the £ = 1,1 = 3/2 phase shift rises quickly from zero to tt, 
and the height of the bump provides a quantitative measure of the overall volume. If the scattering particles 
have intrinsic spin, the denominator in Eq. ([3]) picks up a factor (2Si + 1)(2S , 2 + 1), and the sum over angular 
momenta in the numerator would be expanded, J"^ — > Y] e g(2ff + 1)> where S is the total spin. The phase 
shift would then depend on S in addition to its momentum dependence. Since the inclusion of spin does not 
alter any of the relations derived further below, aside from the addition of spin factors, spin will be suppressed 
throughout the remainder of the manuscript. 

It is easy to derive the Beth-Uhlenbeck formula by considering a scattering center inside a large sphere of 
radius R. The partial wave of outgoing momentum q and angular momentum I would satisfy the boundary 
condition, sm(qR + S() = 0, which gives the constraint qR + Si — nir. Taking the derivative gives dn = d8i/ir, 
thus proving the relation. The relation can be easily extended to inelastic interactions by considering eigenphases 
@. Eq. ([S]) is extraordinarily useful as it also provides the second virial coefficient. Furthermore, since 2d5/de 
is also the time delay in the scattering of a wave packet, it relates the extra time spent near a scattering center 
to the change of the density of states, which can be shown to also satisfy the ergodic theorem [7[ which states 
that in a thermalized system, a particle populates the neighborhood of a scatterer proportional to the local 
density of states. Equation © applies even for relativistic motion, or for energy-dependent interactions. 
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The equivalence between Eq.s (pj and ([3]) in the large volume limit can be explicitly demonstrated by con- 
sidering the large volume limit, where the source function in Eq. (fl]) can be replaced with S(r) «s S(r = 0). In 
that case, 

C(q) = l + 5(r = 0) / d 3 r (|0(q,r)| 2 -l) (4) 

= l + 5(r = 0)47r / r 2 drJ2(2i+ 1)^2 {\^r)\ 2 -\4 0) (q,r)\ 2 }, 
J £ 1 r 

where in the partial wave expansion fa are solutions to the one-dimensional radial Schrodinger equation, with 
an asymptotic form ue{r — > oo) = sin(qr + £tt + 5i(q)). To finish demonstrating the equivalence between the 
Koonin equation and the Beth-Uhlenbeck formula we apply a relation linking the integral of the squared partial 
waves with the derivatives of the phase shifts Q, 

d r{\fa(q,r)\ 2 -\$\q,r)\ 2 } = ^, (5) 
where fag is the partial wave in the absence of the potential. Inserting Eq. ((Sj) into Eq. Q , 

C(q) = l + 5(r = 0)^5>+l A (6) 
q e dq 

which is equivalent to Eq. (|2J) for the case that the source is uniform over a large volume, i.e., S(r = 0) — > 1/V. 

The phase-shift-based expression of Eq. ^ is applicable for a large volume, but for small volumes the Koonin 
form, Eq. (TTJ), is required. The equivalence between the two equations in the large volume limit makes it clear 
that one can identify the change of emission probability arising from final-state interactions with either the 
change in the density of normalized eigenstates, or the change of the amplitude for emitting particles from a 
specific location into the final state. This equivalence emphasizes the role thcrmalization plays in justifying the 
Koonin formula, Eq. |T]), and is also often used to test numerical solutions of the scattered wave function. The 
general assumption has been that any potential that reproduces the phase shifts, can be used in the Koonin 
formula, assuming the range of the potential is much smaller than the size of the source being explored. However, 
as will be shown below, this equivalence becomes invalid when the potential becomes energy dependent. 

The equivalence relation for partial waves, Eq. ((SJ, can be derived from the Schrodinger equation, assuming 
an energy-independentpotential, and using a similar approach as to what is used to derive the basic relations of 
effective range theory gES]. The goal of this brief paper is to show how the equivalence relation is modified for 
the case of energy-dependent potentials. The modification will involve augmenting the squared wave function 
by a simple multiplicative factor that depends on the derivative of the potential with respect to q. Relativistic 
motion of a particle interacting with a vector potential, which has an effectively energy-dependent interaction 
when mapped to the Schrodinger equation, will be considered in detail. After deriving the modification to Eq. 
((SJ) below, the classical limit will be considered, where it will be shown that the same correction factor arises 
from considering the probabilistic enhancement classically, in the presence of an energy-dependent potential. 

The derivation begins by considering the following form for the Schrodinger equation, 

-d 2 + + U(q,r)\ fa(r) = ^/(r). (7) 



Here, U(q, r) would equal 2fiV(r) for the non-relativistic case with no energy-dependent interactions, and fj, 
would represent the reduced mass. We use the term energy- dependent interaction to clarify that the q dependence 
in U(q, r) is a function of the asymptotic kinetic energy and the position. The phrase momentum- dependent 
interaction sometimes refers to exactly such interactions. For instance, in (lll.[l2T| the nuclear optical potential is 
expressed as a function of the magnitude of the local momentum as calculated from classical arguments, which 
means that effectively U(q, r) can be considered energy-dependent. If the potential were expressed as a function 
of r and gradients, i.e., U(V, r), it would not satisfy our definition of being energy-dependent. 

Considering the solution to Eq. for two neighboring values of the asymptotic momentum q and q', 

R dr [-fa(q',r)d 2 fa[q,r)+d 2 fa(q',r)faq,r)] (8) 

rR pR 

+ dr [U(q,r) -U(q',r)](j)'{q',r)(j)(q,r) = (q 2 - q' 2 ) dr (/>' \q ,r)<f>(q,r), 
Jo Jo 
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integration by parts combined with keeping only terms linear in q — q' yields 

[-<P'(q', r)dA{q, r) + d r <f>'(q', r)0(g, r)] r=R = (q- q') J* dr 4> 2 {q, r) Uq - ^U(q, r)) . (9) 

Assuming R is sufficiently large to justify use of the asymptotic form of the wave function, <fii{q, r) ~ sin(qr + 
t-n + Si(q)), and after substituting 0(q, r) — cj)(q', r) + (q — q')d q 4>(q 7 r), 

l( R+ %) = I dr ^ (? ' r) [1 ~ WWWfa r)] • (10) 

After subtracting the same quantity with zero potential (and thus zero phase shift), one finds the generalized 
form of Eq. (O relating dS/dq and the wave function, 

dr {\^(q,r)\ 2 [l-(l/2q)d q U(q,r)}-\^ \q,r)\ 2 } = ~. (11) 

The factor [1 — (l/2q)d q ll(q,r)] is the same as found in [? ] when calculating the local density of states in a 
Green's function approach. The expression can be summed over partial waves to find the analogous expression 
for plane waves, where 0(q, r)| 2 will also be augmented by the factor [1— (l/2q)d q U (q, r)]. However, the extension 
to the outgoing plane-wave case relies on the assumption that U(q, r) is independent of £. This assumption is 
often violated in scattering phenomenology for nuclear physics, where different forms of the scattering potential 
are applied for different £. 

One example where a q dependence for U arises naturally is the Klein-Gordon equation for a vector potential 

V, 

[E - V(r)] 2 <j) = (-V 2 + m 2 )<j). (12) 

After setting q 2 = E 2 — m 2 , the Klein-Gordon equation can be equated to the Schrodinger equation, Eq. ([7]), 
with the potential U, 

U(q, r) = 2EV(r) - V 2 {r), (13) 
which, since d q E = q/E, is momentum dependent so the constraint of Eq. (fTTj) becomes 

R dr [\M<l,r)\ 2 {l-V{r)/E)-\$\q,r)\ 2 ] ■ (14) 

Thus, if one uses |^| 2 (1 — V/E), rather than |^>| 2 , to describe the probability enhancement, it will be consistent 
with the correction to the level density. This result should not be surprising, since the density for a scalar 
field with a gauge field, V — eA , behaves as <j){id t — eA )(j), and for an eigenstate becomes \(/>\ 2 (E — V). This 
explicit appearance of the vector potential into physical quantities related to the field is characteristic of gauge 
invariance, which requires derivatives to be modified, <9 M — > <9 M + ieA^. 

The result is somewhat different when one considers the relative motion of two particles of finite mass, mi and 
m%. Although the following expression neglects retardation effects, the analogy of the Klein-Gordon equation 
for the relative wave function can be written as 

[E - V(r)]<f> = (-V 2 + m 2 )cj) + yj (- V 2 + m 2 )0 . (15) 

Solving for V 2 </> allows one to write an energy-dependent Schrodinger equation for tf>, 

-V 2 cf> = q 2 <P + U(q,r)(f>, (16) 
U(q,r) = q 2 \ [{E V(r)f 2(m 2 + m 2 ) + H ^ 



q 2 = -iE'-2(mi+mi.) 



For the limit that V << E, one can keep only the first order terms in V, and U can be written as 2/ie 2 Jr where 
the effective mass n is 

_ E (m 2 - m 2 ) 2 
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In the non-relativistic limit, E = m\ + m2, and one recovers the usual expression for the effective mass, 
/i = m 1 m 1 /(m 1 + m 2 ), which has no momentum dependence. 

The effects of strong interactions can be described by phase-shifts even when Coulomb interactions are present, 
and as can be seen below, an equivalence relation can be derived linking the derivative of the phase shift to the 
integral of \<f>(q, r)\ 2 — \<fio(q, r)\ 2 , where <f>o includes the Coulomb interaction. Assuming the Coulomb interaction 
is small enough to justify using the approximation C/coui = 2/i(q)e 2 /r, the asymptotic wave function becomes 

4>e(q,r) ~ sin(gr + 7?ln(r) + S e (q)) , (18) 
V = He 2 /q. 

One can repeat the steps used to derive Eq. (fTTj) . but keeping in mind that 4>^ ' was solved with the Coulomb 
interaction. This results in the relation 

dr {\Ml,r)\ 2 [1 - (l/2q)d q U s (q,r) {e 2 /qr)d ql x] - \cf>f(q,r)\ 2 [l - (e 2 /qr)d q ^}} = \^ . (19) 

Thus, even if the short-range potential U s (q, r) has no momentum dependence, the Coulomb interaction modifies 
the equivalence relation by the incorporation of the simple factor [1 — (e 2 / 'qr)d q fj], which has a momentum 
dependence due to the fact that jit depends on q once the motion becomes relativistic, as shown in Eq. (jTTJ) . 

The modification factor [1 — (e 2 /qr)d q p] also appears when one considers the classical analog of the wave 
function. Classically, particles of outgoing momentum q are enhanced by the factor describing the focusing of 
phase space [H], [H| • In these previous studies the role of the squared wave function was shown to be 

|0(q,r)| 2 ^^, (20) 

where q$ is the momentum at the initial separation r. The equivalence works whenever qr ^> 1. Non- 
relativistically, energy conservation implies 

q 2 = q 2 + 2f*j = q 2 o + U(q,r) . (21) 

If one ignores the g-dependence of U, qodqo will equal qdq. Indeed, this independence was assumed when solving 
for the wave function at fixed q, so the classical-quantum equivalence can be stated as 



<l0(q,r)| 2 )^ = Jl-^, (22) 
x q V q z r 

where the (...) denotes an average over angles so that d 3 q — 4irq 2 dq. This will be identical to the squared wave 
function if one solves the Schrodingcr equation for a fixed mass fj,(q) in the limit qr 3> 1. However, the factor 
q a dq /qdq =/= 1 if U depends on q. In that case differentiating the energy-conservation relation gives, 

^L = l-(l/2q)d q U(q,r), (23) 
qdq 

which is the identical factor used to modify the equivalence relation in Eq. (jlip . This emphasizes the physical 
origin of the modification factor. It also suggests that for most applications involving the application of relative 
wave functions for relativistic q, one should incorporate the factor [1 — (l/2q)d q U] into \(j>\ 2 if it is to be interpreted 
as a probability enhancement. 

The Coulomb modification factor above only comes into play for relativistic motion where q is not much 
smaller than the total energy. In the limit of large q, qr is also large and one can justify the classical expression. 
Considering the equal mass case of Eq. (fTF]) . where fi — E/4, the large-g limit of the phase space focusing factor 
becomes 

?r«i-*£--- ( 24 ) 

d a q q z r /ir 

The last term, which arises from the q dependence of /i, is negligible for small q where the motion is non- 
relativistic, but is four times larger than the second term in the relativistic limit. Also of note, both terms fall 
as 1 Jq asymptotically, in contrast to the 1 jq 2 behavior of the non-relativistic expression. 

For most applications, the energy-dependence of the effective mass can be ignored. This is certainly the 
case for low-energy reactions. It is also the case for the vast majority of applications involving two-particle 
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correlations for high-energy collisions, since most analyses focus on correlations at small relative momentum. 
At large relative momentum, corrections for the q dependence of the potential become important, but that 
is also the region where interactions matter less due to the large competing phase space associated with free 
motion. For instance, these effects could be important for modeling tttt interactions in the region of the p 
meson (l6l . [h7| . where the decaying pions are highly relativistic. However, for sources much larger than one 
femtometcr, qr ^> 1, one can model observables using statistical arguments based on derivatives of the phase 
shifts, and avoid applying \(f>\ 2 weights. Another class of observables where these effects might be important are 
analyses of charge balance functions (l8l. [l9|. which involve like-sign subtractions using many pairs, especially 
in a high-energy heavy ion collision. In this case, although the effect of Coulomb interactions at large q is small 
on a pair-by-pair basis, Coulomb effects are magnified by the large number of pairs in a high multiplicity event, 
so that the more complicated structure of Eq. is required. Energy dependent interactions are also applied 
to other forms of reaction theory, e.g., the optical models of [HI, El- Since both the phase shifts and modified 
density are also of interest in many of these applications, the modification factor could also be relevant. In these 
cases the energy dependence is not a relativistic effect, and the characteristic momentum scale for the energy 
dependence is not the rest mass, but instead is determined by some other phenomenological scales. Thus, the 
importance of the modification factor must be assessed on a case-by-case basis. 
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